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Abstract: If scale invariance is exact, unparticles are unlikely to be probed in collid- 
ers since there are stringent constraints from astrophysics and cosmology. However these 
constraints are inapplicable if scale invariance is broken at a scale /i > 1 GeV. The case 
1 GeV < fi < Mz is particularly interesting since it allows unparticles to be probed at and 
below the Z pole. We show that \i can naturally be in this range if only vector unparticles 
exist, and briefly remark on implications for Higgs phenomenology. We then obtain con- 
straints on unparticle parameters from e + e~ — > cross-section and forward-backward 
asymmetry data, and compare with the constraints from mono-photon production and the 
Z hadronic width. 



1. Introduction 



Unparticle physics was introduced in Ref. Q as a low energy effective description of a 
hidden sector with a nontrivial infrared fixed point. This sector is assumed to interact 
with the Standard Model (SM) through the exchange of particles at a high scale M. Below 
M , the interactions are of the form 

(; O^Ouv, (1.1) 



Jyjd uv +d 



SM 



where Cj are dimensionless constants, O l SM is an operator with mass dimension d l SM built 
out of SM fields and Ojjv is an operator with mass dimension duv built out of the hidden 
sector fields. Scale invariance in the hidden sector emerges at an energy scale A < M. In 



the effective theory below A the interactions of Eq. (1.1) take the form 

M d uv+ d hM -i°SMO , (1.2) 

where d is the scaling dimension of the unparticle operator O. 

Unparticle effects might be detectable in missing energy distributions and interference 
with SM amplitudes (l[ |2], ||, [I) . However, if scale invariance is exact, unparticles are unlikely 
to be probed in colliders since there are strong constraints from astrophysics and cosmology 
H U (see Section |2|) . As discussed in Section |3[ these constraints are inapplicable if scale 
invariance is broken at a scale /i > 1 GeV, while constraints from experiments at center-of- 
mass energy yfs > /i remain relevant and resonance-like behavior at \i is expected. Ref. 
has considered collider phenomenology for fi > Mz- Here we consider the constraints on 
unparticle parameters assuming 1 GeV < fj, < Mz which allows unparticles to be probed 
by s channel Z exchange observables. 

For scales of A and M that are experimentally accessible, the Higgs coupling to scalar 
unparticles generally breaks scale invariance at the electroweak scale [||, ||] . Having \i < Mz 
in this case requires somewhat small dimensionless couplings (Section |3.1| ). However, if only 
vector unparticles exist, scale invariance is broken by higher dimensional operators, and [i 
can naturally be below Mz (Section [T^). We also briefly discuss how vector unparticles 
could affect Higgs phenomenology in Section |3.3| . 

Constraints on vector and axial- vector unparticle couplings obtained using e + e _ — > 
forward-backward asymmetry (FBA) and total cross-section data are presented in 
Section |I[ The resonance-like behaviour at /x is taken into account in the analysis. Our 
summary is followed by four appendices covering details of the unparticle contribution 
to the Z hadronic width, the bound from SN 1987A cooling, the vacuum polarization 
correction to the unparticle propagator, and the initial state QED corrections. 

2. Bounds on vector unparticle interactions 

Consider vector unparticles coupling to fermions: 

\dijy— d _ sd uv —d 
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which, following the convention of Refs. [0, ITol, can be written as 



TPT^O^ + T^T^^ , (2-2) 



with 



C ^ = C ^UJ brj • (2 - 3) 

Using the spectral density p{m 2 ) = Ad(m 2 ) d_2 [£]], the propagator is |2], ||, [| 

Here (— q 2 ) d ~ 2 is defined as |<7 2 | d_2 for negative q 2 and |g 2 2 e — * d7r for positive q 2 . Ad is 
chosen following the convention of Ref. Q : 

d (2vr) M r(d-l)r(2(i) ' K ' 

The constant a = 1 if is assumed to be transverse, and a = 2(d — 2) /(d — 1) in conformal 



field theories |TIJ. The value of a does not affect the results of this paper. 

It should be noted that operators of a conformal field theory are subject to lower 
bounds on their scaling dimensions from unitarity, and in particular d > 3 for vector 
operators ]lT| , O]. However, this bound can be violated for a hypothetical scale invariant 
field theory that is not conformally invariant (see e.g. Ref. |R|). We focus on the range 
1 < d < 2 since unparticle effects are relatively suppressed for higher values of d. (Also, 
SM contact interactions induced by messenger exchange at the scale M generally dominate 



over unparticle interference effects for d > 3 [11]. 



A bound on the scale of interactions can be obtained from mono-photon production 
(e + e~ —* 7 + unparticle) at LEP2. The cross-section is given by ||] 

o_ : iJi/r2z7' „ \ A/r2 / 1 „„„2 a 7 ' v / 



8vr 3 M 2 £ 7 s V Mf / 1 - cos 2 6> 7 



where c = \jCy + c Al is the photon energy, and # 7 is the polar angle. Following 
Ref. @, we obtain an upper bound on c using the L3 95% C.L. upper limit a ~ 0.2 pb 
(obtained under the cuts £ 7 > 5 GeV and | cos<9 7 | < 0.97 at ^/s = 207 GeV) @. This 
"mono-photon bound" corresponds to c < 0.026, 0.032 and 0.057 for d = 1.1, 1.5 and 1.9 
respectively. Note that since A < M and unparticle effects can only be probed if y/s < A, 
c > (Mz / yfs)^ 1 is theoretically inaccessible. This implies that the current bound from 
mono-photon production is only relevant for d < 2.6 (see Fig. [l]). 1 

Another process considered in Ref. [H] is Z —* + unparticle, which contributes to 
the Z hadronic width. Here we note that it is important to consider the vertex correction 



1 Mono-Z production is also considered in Refs. [l^]. Similarly to mono-photon production, upper 
bounds on c can be obtained using the L3 limit on Z+ missing energy cross-section , but they are weaker 
than the mono-photon bounds. 
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Figure 1: Upper bounds on c from mono-photon production (solid blue curve), BBN (dashed red 
curve) and SN 1987A (dotted magenta curve). The shaded region corresponds to the theoretically 
forbidden region c > (Mz / \fs) d ~ 1 . 



together with the real emission process, since the two contributions largely cancel each 
other for values of d close to 1 and the former contribution dominates for values of d close 
to 2. As explained in Appendix [A|, the constraint on unparticles from the Z hadronic width 
is also weaker than the mono-photon bound. 

We now compare the mono-photon bound with the constraints on vector unparticles 
from cosmology and astrophysics ||. To preserve the successful predictions of Big Bang 
Nucleosynthesis (BBN), we require the unparticle sector to be colder than SM radiation 
during BBN, so that its energy density is subdominant. For the operator in Eq. fl2.2| ), the 
interaction rate redshifts more slowly than the Hubble parameter H if cZ < 3/2. The 
unparticle sector can then remain cold if it is decoupled throughout BBN, corresponding 
to < Tt for T ~ 1 MeV. For d > 3/2, T^ redshifts faster than TL. In this case we require 
the unparticle sector to decouple before T ~ 1 GeV so that the QCD phase-transition only 
heats up SM radiation ||. The BBN constraint, corresponding to 



( 2 < T*- 2d { T - 1 MeV if d < 3/2 

\M d f l J ~ 10 18 GeV \T~1 GeV if d > 3/2 ' 1 ' ' 

is much more stringent than the mono-photon bound (see Fig. IT]). The SN 1987A constraint 



on unparticle emission p, y, 17 



where the supernova core temperature is taken to be T$n = 30 MeV and ~ 0.01 (see 
Appendix ITJ), is similar in magnitude to the BBN constraint. 
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3. Broken scale invariance 



The BBN and SN 1987A constraints can be evaded provided scale invariance is broken at 
a scale \x sufficiently large compared to the relevant energy scales (~ 1 MeV and ~ Tsn 
respectively). 2 We can model broken scale invariance by removing modes with energy less 
than n in the spectral density, so that || 

p(m 2 ) = A,Mm 2 - f< 2 )(m 2 - M 2 )"- 2 , (3.1) 
[W>U = ^[-tf-^]"(-*. + .^) • (3-2) 

where [— (q 2 — fi 2 )] d ~ 2 is defined as \q 2 —fi 2 \ d ~ 2 for q 2 < and \q 2 — fi 2 \ d ~ 2 e~ ld ^ for q 2 > /j, 2 . 

Due to Boltzmann suppression of the emission, the SN 1987A constraint with scale 
invariance broken at a scale /i corresponds to replacing Cd in Eq. (|2.8|) by (see Appendix 



a / \ d+5/2 

2 m)-d 7r 7/2( d _ 1) 2-d{ TsN ) ■ V'*) 

Assuming c is close to the mono-photon bound, the SN 1987A constraint can be evaded 
provided /j > 1 GeV. Note that other constraints arising from long range forces p0| , 
contributions to the muon and electron anomalous magnetic moments ||, ^] , modifications 



to positronium decay [21|, neutrino decay into unparticles [^], and contributions to low 



energy neutrino-electron scattering amplitudes [23] are also evaded in this case. 



The mono-photon bound is also modified when scale invariance is broken: In Eq. (|2.6|), 
the numerator inside the parentheses is replaced by s — 2y/~sE 1 — fj, 2 , and the end-point for 
E~ is shifted from y/s/2 to (s — /U 2 )/2y / i. However, these modifications do not change the 
cross-section appreciably for fj, < Mz- 

Whether scale invariance is broken or not is relevant for the allowed range of the vector 
unparticle scale dimension d. Consider the decay width from the interaction Eq. ( |2.2j ) 
where an initial fermion with mass mt decays into a massless fermion and the unparticle. 
Following Ref. 24], we obtain 

dT __ A d c 2 E 2 [(2 + a)m 2 -4m f E} 
dE ~ 4vr 2 M| d " 2 (m 2 - 2m f Ef- d 



9{m f - 2E) , (3.4) 



where E is the energy of the final fermion. Integrating over dE, it follows that the total 
decay width diverges for d < 2. This is due to the extra {1/q 2 ) factor associated with the 
vector propagator. However, once scale invariance is broken, values of q 2 < fj? are removed 
from the phase space: 



dT_ A d c 2 E 2 [(2 + a)m 2 - 4m f E] 

dE ~ 4vr 2 M| d " 2 {m 2 - 2m f E)(m 2 f - 2m f E - /i 2 ) 



9(mj - 2m f E - fi 2 ) , (3.5) 



2 Although unparticles are stable if scale invariance is exact, it is not clear if they remain so when scale 
invariance is broken. If they are stable, and if fi is less than the top quark mass, it is not sufficient that 
they decouple at ~ 1 GeV for d > 3/2. Instead, they should remain out of equilibrium at all temperatures 
before BBN, at least up to the reheating temperature |19|. 
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and the total width remains finite and positive for d < 2. 3 

Next, we discuss how scale invariance could be broken such that 1 GeV < \i < Mz, first 
considering the influence of scalar unparticles and then assuming only vector unparticles 
couple to the SM. 

3.1 Scalar unparticles 

As pointed out in Ref. ||, scale invariance is broken by the operator 

A d uv -d 

C2^ji -H^HO, (3.6) 

where H is the SM Higgs doublet, at an energy scale 



/ A \ d uv -2 

^ (£) A " 



l/(4-d) 

(3.7) 



where v = 174 GeV. Having an experimentally accessible conformal window ji <C A ~ v 
requires C2 < 1 H. Assuming fi < Mz, another upper bound on \i and C2 can be obtained 
from the threshold correction to the fine structure constant @. If the operator 

\d uv -d 

c ±l^r FPF ^ < 3 - 8 ) 

exists, the value of a~ 1 (Mz) remains within the current uncertainty for 

f M \ duv/d f 1 \ 1/d A , , 



Eq. ( |3.9D provides an upper bound on /j,, whereas Eqs. ( |2.8| , p^ ) provide a lower 
bound. There can be a scale invariant window below Mz between these two bounds 
without violating any of the other constraints discussed above. As a specific example we 
take duv = 2 or 3, Cy = Ca = 1 and A = v. Setting M equal to the mono-photon 



bound using Eqs. (|2.3j , pTq ), we calculate the range of C2, C4 and \x that satisfies the other 
constraints. As shown in Table |] and Fig. ^, there is an allowed range of [i below Mz, 
provided the scalar unparticle operators couple somewhat weakly (C4, C2 < 0.1) compared 
to vector operators {Cy = Ca = 1)- 

3.2 Vector unparticles 

Even if only vector unparticles exist, scale invariance can still be broken if the Higgs 
couples to higher-dimensional operators such as O^O^. Furthermore, due to the higher 
dimensionality, the scale fi is naturally suppressed compared to the electroweak scale. 4 
Consider the operators 

\2d uv ~d* \2d uv -d, 

—g-^fiHOPOr + ——.f^f^o^ (3.10) 

3 For scalar unparticles and d < 1, the divergence pointed out in Ref. |l| remains whether fi — or not. 
4 See Ref. pa] for a similar scenario with charged scalar unparticles. 
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djjv — 2 duv — 3 



d 


M (GeV) 


fi (GeV) 




C 2 


M (GeV) 


H (GeV) 




c 2 


1.1 


8800 


1.5-M Z 


9 x 10" 7 - 


-0.15 


1200 


1.5-35 


7 x 10" 6 


-0.07 


1.5 


5600 


1.3-74 


4 x 10" 6 - 


-0.12 


990 


1.3-31 


2 x 10" 5 


-0.08 


1.9 


2400 


1.1-49 


2 x 10" 5 - 


-0.07 


650 


1.1-25 


8 x 10" 5 


-0.06 



Table 1: The allowed range of ^ (below Mz) and C2 for M at the mono-photon bound, assuming 
C v = C A = 1 (see Eq. (pT|)), A = u and C 4 = C 2 (see Eqs. |Dj, pD). 
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Figure 2: The shaded region in the C4-C2 plane is allowed for M at the mono-photon bound 
(Cy = CU = 1, A = v). The lower horizontal line corresponds to the SN 1987A constraint, the 
middle horizontal line corresponds to /i = Mz, and the upper line to /i — v. The diagonal line is 
the constraint from a(Mz). 



where we have set C2 = C4 = 1, and the scale dimension of O^O^ = < 2<i. Eq. (| 
and Eq. (|3.9[) are modified as follows: 



(£) 



l/(4-d») 



M 

T 



2d uv /d* 



10" 4 - 5 / d *A. 



(3.11) 
(3.12) 



As shown in Fig. 0, ^ can easily lie in the allowed range. 
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Figure 3: The shaded regions in the d-fj, plane are allowed for M at the mono-photon bound 
(A = v, C v 



C, 



1). The solid blue curve is /i as given by Eq. (3.11). The dashed red 



curve corresponds to the SN 1987A constraint, the dotted-dashed magenta curve corresponds to 
the constraint from a(Mz), and the horizontal dotted line corresponds to fi = Mz- 



3.3 Implications for Higgs phenomenology 

The effects of scalar unparticles on Higgs phenomenology have been considered in Refs. 
1 26, 27 1. 5 For scalar unparticles the same operator H* HO is responsible for breaking scale 
invariance and Higgs-unparticle mixing to lowest order, and thus the effects are suppressed 
for [i <C Mz- To be more explicit, the mixing between the SM Higgs boson h and the 
unparticle is induced by the interaction term ((i A ~ d /v)Oh. Considering the effective Higgs 
coupling (l/v)C~ / r y hF fll ,F tlL ' as an example, the contribution from the above interaction and 
Eq. is given by |26] 

c^o^c^(JLrm d (^r° r . <3. 13 ) 



77V '" " ' ' " 2sin(f7r \m h J \AJ \M J 

Provided fi <C Mz, this is small compared to the SM effective coupling C 77 ~ 10 -3 . 
The Higgs partial decay width to fermions induced by the operator ip^^D^ipO is similarly 
suppressed. 

On the other hand, for vector unparticles scale invariance is broken by H^HO^O^ 
whereas mixing is (also) induced by 



M d w 



^H^D^HO 1 * . (3.14) 



5 See also Ref. [E8| for supersymmetric unparticle effects. 
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Using Eq. (|2.1| ) and Eq. ( [3.14 ), the effective fermionic operator is 

-LrfD^H^, (3.15) 



where 



A 2 



1 e- id *A d (J^\ 2d ~ 2 ( A ^ Ww-V 



A 2 ff s 2 sin dir \ A / \M 



(3.16) 



Contributions of this operator to Higgs production at a linear collider have been considered 
in Ref. [E9]. The main effect is interference with the SM Higgs-strahlung (HZ) cross-section, 
which can be substantial in the e + e~ — ► h[i + n~, Jit + t~, hqq channels for M close to the 
mono-photon bound. 

It is also interesting to note that the operator HO^O^ induces a partial decay width 
T(h — > O^O^) ~ / u 8-M */(t> 2 m^ -M *). Although this is typically small, it becomes of order 
m\/v 2 in the limit — ► 4 (i.e. d* — ► 2d and d — > 2). With decays of 0^0^ suppressed, 
the Higgs would then decay invisibly. 



4. Muon pair production bounds on vector unparticles 



We have already obtained a collider bound using Eq. (|2.6|). Other bounds can be obtained 
using the ratio Rjj = <r(with unparticles) ja (without unparticles) as well as the FBA (de- 
fined in Appendix O) for e + e~ —* fi + fi~, and by combining measurements at and away 



from the Z pole. As shown in Ref. y] and discussed further in Ref. [ 10 1 , vector cou- 
plings of unparticles will mainly affect Ru away from the Z pole, and FBA at the Z pole. 
Axial-vector couplings have the opposite behaviour. 

Due to the resonance-like behavior at fi (referred to as "un-resonance" 0]), measure- 
ments at energies around fi would be particularly sensitive to unparticle effects. Thus the 
bounds on cy,A (defined in Eq. ( |2.2| )) for a given value of d will also depend on fx. As 
an example we plot FBA and Ru for d = 1.1 in Fig. ^. Taking ca = 0.026 and cy = 0, 
FBA= -7.2% for ^Ts = 34.8 GeV if \i = 30 GeV, to be compared with -8.3% if fi < 30 
GeV, and —8.9% for SM. Taking into account the measurement FBA= — 10. 4 ± 1.3 ±0.5% 
at the same center-of-mass energy p0| , it is clear that the bound on ca for fj, = 30 GeV 
will be more stringent compared to the bound for \i <C 30 GeV (see Fig. ||). 

It should be noted that for the propagator in Eq. ( |3.2| ), the area under the un-resonance 
diverges for d < 1.5. However, it is likely that once scale invariance is broken, particle-like 
modes will appear in the spectral density ||]. For example, vacuum polarization correction 
from fermion loops will modify Eq. ( |3,2j ) as follows: 

1 i 

(4.1) 



(q 2 - /i 2 ) 2 - rf (q 2 - fi 2 ) 2 ~ d - n(q 2 ) 

It can therefore be expected that the unparticle will become unstable, and the area under 
the un-resonance will depend on the decay width. 

We have performed a x 2 analysis of LEPl-Aleph, KEK- Venus and PETRA-MarkJ 
e + e~ — > cross-section and FBA data p0, En]. The simulation includes the vacuum 
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vector coupling 



axial-vector coupling 
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Figure 4: FBA and Ru for e + e _ — > with d = 1.1. Solid green curves: SM; dashed blue 

curves: unparticlcs with /i — 0, dotted red curves: unparticles with /i = 30 GeV. (cy,.A = 0.026 
correspond to the mono-photon bound of section ^[) 



polarization correction from fermion loops to the unparticle propagator (see Appendix |C|) 
and uses a fixed Z decay width Tz = 2.41 GeV which is the SM best-fit value for the data. 
Initial-state QED corrections are also included (see Appendix [D]). 

The allowed regions in the cy-CA plane for different values of d and \i are shown in 
Fig. [5]. The best-fit parameters and x 2 values are listed in Table |2], and fits to FBA data 
with and without unparticles are displayed in Fig. |6[ For values of d close to 1 where 
fermion- unparticle couplings are less suppressed by M^~ d , constraints on cy and ca are 
more stringent and the dependence on \i is more significant. The mono-photon bound 
discussed in Section § is stronger than the muon pair production bound for d > 1.3. 6 

Finally, it is worth emphasizing that the spin and scale dimension of the exchanged 
unparticle can be probed by analyzing the scattering angle and energy distributions of 
differential cross-sections in a linear collider, for both real emission and virtual exchange 
processes (2|, |||, [Tq] . Furthermore, for polarized beams the azimuthal dependence of the 
final state fermion can provide an independent measure of the scaling dimension for spin-1 
unparticle exchange 

6 Recently, it was noted that processes mediated by unparticle self-interactions lead to multi-body final 
states which could be the most promising modes for unparticle discovery at colliders |}2| . However, details 
of the hidden sector are required to make predictions. 
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Figure 5: Allowed regions in the cy-CA plane from a \ 2 analysis of e + e" — > cross-section 
and FBA data. The contours represent the la, 2a and 3a regions. We only show results in the 
first quadrant since the dependence of the cross-section and FBA on the relative sign of cy and Ca 
is too weak to be visible. 
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fx = 1.5 GeV n = 30 GeV fx = 70 GeV 
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0.0081 


0.020 


154.3 


1.5 


0.0089 


0.093 
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0.0059 


0.081 


153.6 


0.0016 


0.036 


154.2 


1.7 


0.083 


0.13 


153.1 


0.12 


0.12 


152.0 


0.0071 


0.081 


153.9 


1.9 


0.11 


0.11 


152.6 


0.12 


0.11 


152.1 


0.085 


0.11 


153.4 



Table 2: Best-fit parameters and \ 2 values from an analysis of LEPl-Aleph, KEK- Venus and 
PETRA-MarkJ e + e~ — > cross-section and FBA data. The dataset is comprised of 55 FBA 

data points and 54 cross-section data points. The x 2 value for SM is obtained from a scan on 4 SM 
parameters Mz, ^z, e and sm6\y. The \ 2 values for SM + unparticles are obtained from a scan 
on 2 unparticle parameters cy and ca for different fixed values of d and fi, with SM parameters 
fixed to their SM best-fit values. 
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Figure 6: Fits to e + e~ — > [i + [iT FBA data, from a scan over cy and ca for different fixed d and 
ji values. The solid curve (red) is the SM fit and the dashed curve (blue) is the best-fit curve for 
d = 1.9 and /j, = 1.5 GeV, /i = 30 GeV or /i = 70 GeV. For d < 1.5 the curves with unparticle 
contribution are almost indistinguishable from the SM curve. Note that the unparticle FBA curve 
does not exhibit divergent behavior at \i as the vacuum polarization introduces a finite decay width 
and stabilizes the unparticle propagator. 
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5. Summary 

For exact scale invariance, astrophysical and cosmological constraints are in gross conflict 
with the possibility of probing unparticles in colliders. We showed that for vector un- 
particles collider constraints become relevant only if scale invariance is broken at a scale 
fj, > 1 GeV. Breaking the scale invariance also affects collider expectations by giving rise 
to a resonance-like behaviour. On the other hand, unparticle effects cannot be observed 
at energies below the scale [i. We focused on the case 1 GeV < fi < Mz which allows 
unparticle effects to show up in Z exchange observables, and gave demonstrations of how 
this can be realized through unparticle-Higgs couplings. 

Simple bounds on vector unparticles have been obtained using effective contact inter- 
actions in Refs. Q [|. Here we have made a more detailed analysis using e + e~ — > /U + /i~ 
cross-section and forward-backward asymmetry data both at the Z pole and away from 
it, also taking into account the resonance-like behaviour associated with broken scale in- 
variance. We found that unparticle parameters are severely constrained for values of scale 
dimension d close to 1. For d > 1.3, constraints from mono-photon production are more 
stringent compared to constraints from muon pair production. 
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A. Unparticle contribution to the Z hadronic width 



Ref. 1 34] studied the real and virtual massive vector boson contribution to the Z hadronic 
width Rz- To calculate the constraint on unparticles, we write the unparticle operator in 
terms of deconstructed particle fields [35|: M = YljFj^ji where the field has mass 
Mf = j A 2 and 



F 2 = ^A 2 (M, 2 ) d - 2 
J 2ir J 



(A.l) 



In the limit A —* 0, the contribution to Rz is obtained by integrating the contribution 
from a vector boson with mass m over 5 = m 2 /M^: 



AR Z = A d c 2 v 
Rz ~ 16vr 3 



5 d - 2 F l {8)d8 + 



F 2 (5)d5 



(A.2) 



where [34| 



F 1 (S) = {l + 5f [3 In 6 + (In 5) 2 ] + 5(1 - 5 2 ) - 25ln5 



-2(1 + 5) 2 



ln(l + <5) In5 + Li 2 



l + <5 



Li 2 



l + <5 



(A.3) 
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F- 2 (8) - -2{ 7 - + 5+(5+^)ln5 



Li 2 T-^ + 7: In 



1 , of 8 \ 7T 



1 + 5/ 2 V 1 + 5/ 6 



^2 



(A.4) 



Li2(x) = — Jj^ diln(l — is the Spence function, and uniform coupling for quarks is 
assumed. Note that the upper limit 1 of the 5 integration is kinematic for the real emission, 
and the upper limit becomes oo for the virtual correction. 

Evaluating the integrals, we obtain ARz/Rz — O.Olcy, corresponding to a bound 
cy < 0.3 since ARz/Rz = Aa s /n ~ 0.001. Including the axial-vector coupling is straight- 
forward and leads to c < 0.3/\/2. 

B. The bound from SN 1987A cooling 

As discussed in Refs. §, |, |||, ||], SN 1987A energy-loss arguments provide very 
restrictive constraints on unparticle couplings. In this section we discuss the constraint 



from pair annihilation of neutrinos and obtain the prefactor Cd in Eqs. ( [2.8] , 3.3) following 
the method in Refs. pL M. 7 



The observed duration of SN 1987A neutrino burst puts a constraint on the supernova 
volume emissivity [[37]] 

Q < 3 x 10 33 erg cm' 3 s' 1 , (B.l) 
where the supernova core temperature is taken to be T$n = 30 MeV. This corresponds to 

Q < 4 x lO- 22 ^ . (B.2) 

As in Appendix |A|, we write the unparticle operator in terms of deconstructed particle 
fields. The cross-section for neutrino pair annihilation to \j is 

/ \ 2 

A d A 2 (Mf) d - 2 5(s-Mf). (B.3) 



0i 



mJ" 1 



The supernova volume emissivity is found by thermally averaging over the Fermi-Dirac 
distribution (see e.g. Ref. 



Q ' - 1 ^i^i+i / j^k^hi {Ei + Eh)2 " Ti • (B - 4) 

where we ignored chemical potentials (see Ref. [17]]), and s = 2E\E2(1 — cos#). The total 
emissivity is obtained as 8 



A* J ^^-^{Mt 1 ) SN ' 



7 The constraint from pair annihilation (for exact scale invariance) is discussed in Ref. The con- 

straint from nucleon bremsstrahlung is similar in magnitude ^| [l7|] . 
8 See Ref. fed] for similar calculations with tensor unparticles. 
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where 

2 2d-3 Ad roc (x 1 X2) d (x 1 +X 2 ) 



Cd = jt^ / dxxd^ ^ ^ \ x *' ~ 0.01 . (B.6) 

vr 4 d y (e Xl + l)(e X2 + 1) 

We now repeat the calculation for non-zero \x. By matching to the spectral density in Eq. 
Q3.1| ), we have 

F j = ^-A 2 (M 2 - ff-H{M] - fj, 2 ) , (B.7) 



2 



C 



— AdA (Mf - ^)^0{Mf - ^)5{s - Mj) . (Bi 
M z J 



Using Eqs. ([Bj, |B^8|), we obtain Eq. ([BJ) with 



2 d ~ 3 Ad f°° f°° /" 1 -A t2 /(2r| Jv a:iX2) 

Cd = t — / dxi / dx 2 / d(cos6>) 

Jo 7 M 2/(4T2xi) 7-1 



(xix 2 ) d (xi + x 2 )(l - cos 0) (l - cos# - 2T / 2 r V " 



(e Xl + l)(e X2 + 1) 

The approximation Eq. ( |3.3[) is obtained from Eq. ( |B.9[) assuming /u 3> T$n- 



C. Vacuum polarization correction 

To lowest order, H(q 2 ) in Eq. (f4.1|) is given as follows: 



(B.9) 



n = n LL + u LR + n flL + , (c.i) 

n - - =» - - 2 ^S^!? x ^ US-** ) • (&2) 

n ""»< = - 4 si|^ i 1 di - 1)92 - H ^ ( m »-^i-x tf ) ' <a3) 

where C£ = cy — c^, cr = cy + ca and m/ is the mass of the fermion in the loop. II(g 2 ) 
is complex for the s channel with q 2 > 4m 2 , and the imaginary part will stabilize the 
propagator when the real part coincides with the pole. We assume a universal coupling 
between the unparticle and different fermions that include charged leptons, neutrinos and 
quarks. A numerical example for H(q 2 ) that is summed over the fermions is shown in 
Fig. ft 



D. Initial state QED corrections 

Initial state QED corrections significantly affect the cross-section and FBA around \x (see 
Fig. H). Since the corrections to the SM cross-section osm are removed from the KEK- 
Venus and PETRA-MarkJ data, we only consider the corrections to the unparticle exchange 
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Figure 7: n(g 2 ) from charged lepton, neutrino and quark loops, assuming fermion couplings 
Cy = ca = 0.05. 
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Figure 8: FBA and Ru for e + e~ — > /i + /.t~, d = 1.1 and /i = 30 GeV. Solid blue curves: with initial 
state QED corrections to the unparticle exchange term as well as the interference terms between 
7, Z and the unparticle. Red dashed curves: without initial state QED corrections. (cv,a — 0.026 
correspond to the mono-photon bound of section |[) 



term o\j and the interference terms between 7, Z and the unparticle. The corrected 
cross section is obtained by convoluting the relevant terms with a radiator function H(x): 



o-(fl) = (J SM {s) + / M dxH(x) (au[s(l - x)\ + a int [s(l - x)\) , (D.l) 
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where |J9| 



H(x) = (3x l3 - 1 5 v + 5 h 



(D.2) 



with 



5 h 



-(L-l), L = log4 
~(L-l)(z-2) + ... . 

7T 



1 + 



a ( 3 



vr V 2 



-L + 



7T 



2 +... , 



(D.3) 



The LEPl-Aleph data are fitted with full QED corrections, since the corrections to cjsm 
are not removed. 

The corrected FBA for KEK- Venus and PETRA-MarkJ data is obtained in a similar 
manner IfiOl: 



FBA(s) 



1 



a(s) 



dz 



4z 



H(z){o5 B (z8) + o££(zs)) 



FBi 



Aml/s (1 + zf 



where 



a 



FB 



6>it/2 



do_ 
'dJ2 



da 



(D.4) 



(D.5) 



H(z) = H(l -z)+ (|-) L 2 



(1 ~ zf 
2z 



;i + z)io g (z) + 2(i-z) 



+ 



2 ' 1 

arctan — = — arctan 



(D.6) 



Again, the LEPl-Aleph data are fitted with full QED corrections. 

At the energy scale Mz with the scaling breaking parameter fj, > 1 GeV, unparticle 
bremsstrahlung is not effective and thus not included. 
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